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Abstract 

We have explored a simple microscopic model to simulate a thermally acti- 
vated rate process where the associated bath which comprises a set of relaxing 
modes is not in an equilibrium state. The model captures some of the essential 
features of non-Markovian Langevin dynamics with a fluctuating barrier. Mak- 
ing use of the Fokker-Planck description we calculate the barrier dynamics in 
the steady state and non-stationary regimes. The Kramers-Grote-Hynes reac- 
tive frequency has been computed in closed form in the steady state to illustrate 
the strong dependence of the dynamic coupling of the system with the relax- 
ing modes. The influence of nonequilibrium excitation of the bath modes and 
its relaxation on the kinetics of activation of the system mode is demonstrated. 
We derive the dressed time-dependent Kramers rate in the nonstationary regime 
in closed analytical form which exhibits strong non-exponential relaxation ki- 
netics of the reaction co-ordinate. The feature can be identified as a typical 
non-Markovian dynamical effect. 
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I. Introduction 



More than half a century ago Kramers 1 considered the problem of activated rate 
processes by using a model Brownian particle trapped in a one dimensional well which is 
separated by a barrier of finite height from a deeper well. The particle was supposed to 
be immersed in a medium such that the medium exerts a frictional force on the particle 
but at the same time thermally activate it so that the particle may gain enough energy 
to cross the barrier. Over several decades the model has been the standard paradigm 
in many areas of physics and chemistry 2 . The Kramers problem was to find the rate 
of escape from the well to the barrier. The motion of the particle is governed by the 
following phenomenological Langevin equation, 

1 dV(x) 1 . , 

x = ±±- lx + - F (t) , l 

m ox m 

where x is the coordinate of the particle of mass m moving in a potential V(x). 7 and 
F(t) are the damping rate and the Gaussian stationary random force provided by the 
thermal bath respectively. The properties of noise can be summarized by the following 
two relations, 

(F(t)) = , (F(O)F(f)) = 2 im KT5{t) . (2) 

The Langevin equation (1) is equivalent to the Fokker-Planck equation for proba- 
bility distribution p = p(x,v,t) [also known as Kramers equation], 



dp 1 dV(x) dp dp 



- v— + 7 



dt m dx dv dx 



KT d 2 p d 



m dv 2 dv 



(3) 



Kramers 1 obtained the steady state escape rate k in the limiting cases of high and 
low damping rates in the following form, 

, ^expf— |y 7 — > 00 

k={ 2 " 7 L KTl , (4) 



where u Q and u b are the frequencies associated with the curvature of the potential at 
the bottom of the well and at the barrier top, respectively. Ef, refers to the depth of 
the well. Kramers has also derived an expression for 'intermediate' value of 7 : 



k "° 



2-KUJb 



7 
2 



exp(-E b /KT) 



For non-Markovian random processes where one takes into account of the short 
internal time scales of the system compared to that of the thermal bath, the Langevin 
equation(l) gets replaced by its non-Markovian counterpart 3,4 , sometimes called the 
generalized Langevin equation (GLE); 

1 dV(x) 



x = 



f drZ(t-r)x(r) + —R(t) , (5) 
Jo m 



m dx 

where R(t) is Gaussian but non-Markovian such that 

(R(t)) = 0, (R(0)R(t)) = Z(t)mKT . (6) 
The memory function Z(t) is expressed in terms of Fourier-Laplace components 

Z n {u) = / dtZ(t)e- m ^ (7) 

Jo 

with Zq{uj) = 7 

Based on equation (5) Adelman 5 obtained the generalized Fokker-Planck equation 
for a Brownian oscillator with a parabolic potential as given by ; 

dt UJbX dv V dx ^ dv V ^ ^ m dv 2 m \ou 2 J dvdx 
where 7 = 7(0 and u% = ou 2 (t) are now functions of time [although bounded , they 
may not always provide long time limits] which play a decisive role in the calculation 
of non-Markovian Kramers rate. 

Various workers have made use of generalized Langevin equation to treat the differ- 
ent aspects of the escape problem in the non-Markovian regime. For example, Grote 



and Hynes 4 considered the average motion of the particle in the vicinity of the barrier 
governed by GLE and found that on the average the particle is slowed down by fric- 
tion and defining a reactive frequency A r they showed that the average motion goes as 
exp(±A r t). The analysis of Hanggi and Mojtabai 6 on the other hand is based on the 
generalized Fokker-Planck equation of Adelman with a parabolic potential in the high 
friction limit. The generalized FP approach has also been adopted by Carmeli and 
Nitzan 7 to derive the expression for the steady-state escape rate in the high and low 
friction limit in the Markovian as well as non-Markovian regimes. A comprehensive 
overview has been given in Ref.(2). 

While the early post-Kramers development as summarized above is largely phe- 
nomenological, an interesting advancement in the theory of activated rate processes 
was made when the generalized Langevin equation was realized in terms of a micro- 
scopic model which comprises a system coupled linearly to a discrete set of harmonic 
oscillators. Using the properties of the bath and a normal mode analysis it was shown 8 
that the reactive frequency A r defined by Grote and Hynes 4 for the average motion 
across the barrier is actually a renormalised effective barrier frequency. 

The object of the present paper is twofold : First is to consider a simple variant of the 
system-heat bath model 9 ' 10,11 to simulate the activated rate processes, where the asso- 
ciated bath is in a nonequilibrium state. The model incorporates some of the essential 
features of Langevin dynamics with a fluctuating barrier which had been heuristically 
and phenomenologically proposed earlier in several occasions. 10,13-17 While the major- 
ity of the treatments of the phenomenological fluctuating barrier rest on the reduction 
of the equations to overdamped limit 5,10,14 , thus restricting the validity of the solutions 
in the large time limit, we take full account of the inertial terms in our calculation 
of barrier dynamics and probability distribution function both in the long time and 
in the short time nonstationary regimes. The Fokker-Planck description allows us to 
calculate Kramers-Grote-Hynes reactive frequency pertaining to these situations for 
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non-Markovian dynamics in closed form. Second, since the theories of activated pro- 
cesses traditionally deal with stationary bath, the nonstationary activated processes 
has remained largely overlooked so far. We specifically address this issue and examine 
the influence of initial excitation and subsequent relaxation of a bath modes on the 
activation of the reaction co-ordinate. We show that relaxation of the nonequilibrium 
bath modes may result in strong non-exponential kinetics and a nonstationary Kramers 
rate. The physical situation that has been addressed is the following : 

We consider that at t — 0_, the time just before the system and the bath is subjected 
to an external excitation, the system is appropriately thermalized. At t — 0, the 
excitation is switched on and the bath is thrown into a nonstationary state which 
behaves as a nonequilibrium reservoir. We follow the stochastic dynamics of the system 
mode after t > 0. The important separation of the time scales of the fluctuations of 
the nonequilibrium bath and the thermal bath (to which it relaxes) is that the former 
effectively remains stationary on the fast correlation of the thermal noise. 

The outline of the paper is as follows; Following Ref. [10] we discuss in Sec. II a 
microscopic model to simulate an activated rate process where the system in question 
is not initially thermalized. Appropriate elimination of reservoir degrees of freedom 
leads to a nonlinear non-Markovian Langevin equation which governs the dynamics of 
a particle with a fluctuating barrier, stochasticity being contributed by both (additive) 
thermal noise and a slower (multiplicative) noisy relaxing nonequilibrium modes. The 
Fokker-Planck description is provided in Sec. III. The standard Markovian description 
and the generalized FP equation of Adelman's form can be recovered in the appropri- 
ate limits. In Sec. IV we derive the expression for Kramers rate of barrier crossing in 
the non-Markovian but steady state regime and show that the Kramers-Grote-Hynes 
"reactive frequency" can be explicitly realized in this model in closed form. Sec.V is 
devoted to nonstationary aspect. We solve the time-dependent FP equation for non- 
stationary probability density and calculate the corresponding current. An expression 
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for Kramers rate in the nonstationary regime in closed analytical form is derived. The 
paper is concluded in Sec. VI. 



II. The model and the Langevin equation 



We consider a model consisting of a system mode coupled to a set of relaxing modes 
considered as a semi-infinite dimensional system ({g^j-subsystem) which effectively 
constitutes a nonequilibrium bath. This, in turn, is in contact with a thermally equili- 
brated reservoir. Both the reservoirs are composed of two sets of harmonic oscillators 
characterized by the frequency sets {uik} an d {Qj} for the nonequilibrium and the equi- 
librium bath, respectively. The system-reservoir combination evolves under the total 
Hamiltonian 



the first two terms on the right hand side describe the system mode. The Hamiltonian 
for the thermal and nonequilibrium baths are described by the sets {Qj, Pj} and {qj, pj} 
for coordinates and momenta, respectively. The coupling terms containing Kj refers to 
the usual system-thermal bath linear coupling. The last two terms indicate the coupling 
of the nonequilibrium bath to the system and the thermal bath modes, respectively. 
Since in the present problem, H is considered to be classical and temperature, T high 
for the thermally activated problem we note that quantum effects do not play any 
significant role. Hamiltonian (9) is a simpler variant of that treated in Ref.[10]. For 
simplicity we take m — 1 in (9) and for rest of the treatment. As shown in Ref. [10] 
the model (9) captures the essential features of fluctuating barrier dynamics. We recall 
the relevant aspect in the following discussions. 




+ v(x) + \ + njoj) + \ + uihl) 



~ x H K iQj ~ 9{x) Qk ~ Yl a jkQkQj 



(9) 



j k j,k 
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Eliminating the equilibrium reservoir variables {Qj,Pj} in an appropriate way 9,10 
one may show that the nonequilibrium bath modes obey the following equations of 
motion, 

9* + 7<7fc + ^kQk = g(x) +rj k (t) . (10) 

This takes into account of the average dissipation (7) of the nonequilibrium reservoir 
modes q k due to its coupling to thermal reservoir which induces fluctuations r)k(t) char- 
acterized by (r)k(t)) = and the usual fluctuation-dissipation theorem (r)k(t)r)k(0)) = 
2 r yKTS(t). We mention here that moving from Eq.(9) to (10) generate cross terms of 
the form J2jlkjQj, which are neglected for j 7^ k. 

Proceeding similarly to eliminate the thermal reservoir variables from the equations 
of motion of the system mode one obtains 

x + 7e q i + V'(x) = £ cq (t) + g'(x) Qk , (11) 

k 

where 7 eq refers to the dissipation coefficient of the system mode due to its direct 
coupling to the thermal bath providing fluctuations £ e q(^)- Here we have 

(&q(*)> = and (U(*)U(0)> = ^KT5(t) . 

Now making use of the formal solutions of Eq.(10) 10 which takes into account of the 
relaxation of the nonequilibrium modes and integrating over the nonequilibrium modes 
with a Debye type frequency distribution of the form, 

!?>uj 2 /2ul , for I a; I < u c 
, for I a; I > uj c 

where u c is the high frequency Debye cut-off, one finally arrives at the following 
Langevin equation of motion for the system mode, 

x + r(x)x + V\x) = £ eq {t) + g'{x)£ neq (t) . (12) 
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Here T(x) is a system coordinate dependent dissipation constant composed of 7 eq 
and 7 n eq as follows, 

T(x) = lcq + lncci [g'(x)} 2 . (13) 

£ neq refers to the fluctuations of the nonequilibrium bath modes which effectively cause 
a damping of the system mode by an amount 7 neq [(/(a;)] 2 . 

Eq.(12) also includes the modification of the potential V(x) in which the particle 
moves as 

V(X) = V(X) - ^ Tneq 9 2 {X) . (14) 
71 

Eq.(12) thus describes the effective dynamics of a particle in a modified barrier, 
where the metastability of the well originates from the dynamic coupling g(x) of the 
system mode with the nonequilibrium bath modes. It is necessary to stress here that 
g(x), in general, is nonlinear. This nonlinearity has two immediate consequences. First, 
by virtue of the term V'(x) in Eq.(12) it gives rise to a fluctuating barrier. Second, 
the term g'(x)^ neq (t) imparts a multiplicative noise term in Eq.(12) in addition to the 
usual additive noise term £ e q(£)- We point out here that the problem of diffusion over 
a fluctuating barrier 13-1 of similar nature has been addressed earlier by a number of 
workers from the phenomenological point of view. For example, Stein et.al. 14 have 
calculated the decay of probability from the metastable state in the white noise limit 
and also for short finite correlation times for the fluctuating part of the potential. 
Riemann and Elston 15 have calculated an asymptotic rate formula when the particle 
is subjected to both dichotomous and thermal noise. 

The treatment followed in the aforesaid cases concerns overdamped situation and, in 
general, the validity is restricted to long time limit. In the present problem, however, 
we look at the stochastic process right from the moment the nonequilibrium excitation 
(followed by the relaxation) sets in. We are therefore forced to take into consideration 
of the inertial term in Eq.(12) on its usual footing. 



8 



We now turn to the another aspect of the problem. In order to define the prob- 
lem described by Eq.(12) completely, it is further necessary to state the properties of 
fluctuations of the nonequilibrium bath £ nC q(£)- We have first for Gaussian noise 

(W*)> = • 

Also the essential properties of £ ne q(£) explicitly depend on the nonequilibrium state of 
the intermediate oscillator modes {q^} through U(u, t), the energy density distribution 
function at time t in terms of the following fluctuation-dissipation relation 10 for the 
nonequilibrium bath, 

U(u,t) = / dr (U q (m ncq (t + r)) e^ T 

47neq J-oo 



2 



U(u,0)-±KT 



(15) 



U(u,0) — \KT is a measure of departure of energy density from thermal average 
at t — 0. The exponential term implies that this deviation due to the initial ex- 
citation decays asymptotically to zero as t — > oo, so that one recovers the usual 
fluctuation-dissipation relation for the thermal bath. With the above specification 
of correlation function of £ neq Eq.(15) thus attributes the nonstationary character of 
the {(^-subsystem. 

In passing, we stress that the above derivation 10 is based on the assumption that 
£ neq is effectively stationary on the fast correlation of the thermal modes. This is 
a necessary requirement for the systematic separation of time scales involved in the 
dynamics. We point out that the effective dynamics sets no choice on any special form 
of coupling g(x) between the system mode and the relaxing mode and as such this may 
be of arbitrary nonsingular type for our problem we have considered here. 
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III. The generalized Fokker-Planck description 



Eq.(12) is the required Langevin equation for the particle moving in a modified po- 
tential V(x) [Eq.(14)] and damped by a coordinate-dependent friction T(x) [Eq.(13)] 
due to its linear coupling to a thermal bath and nonlinear coupling to the {qk}- 
subsystem characterized by fluctuations £ ne q(£)- Before proceeding further a few per- 
tinent points are to be noted to stress some distinct and important aspects of the 
model. 

First, depending on the system- {g fc }-subsystem coupling g(x) both the modified 
potential V(x) as well as T(x) are, in general, nonlinear. So the stochastic differential 
equation (12) is nonlinear. Again, the stochasticity in Eq.(12) is composed of two 
parts : £ eq (£) is an additive noise due to thermal bath while £, n eq9'(x) is a multiplicative 
contribution due to nonlinear coupling to {gfc}-subsystem. It is thus important to note 
that the presence of multiplicative noise and a fluctuating barrier are associated with 
nonlinearity in g(x). 

Second, the Langevin equation (12) is non-Markovian. The origin of this non- 
Markovian nature lies in the decaying term in Eq.(15) where the decay explicitly ex- 
presses the initial nonequilibrium nature of the {g/cj-subsystem following the sudden 
excitation at t — 0. This non-Markovian feature is thus not to be confused with that 
arises due to the usual frequency dependence of the dissipation constant. 

Third, although the modification of V(x) is due to the specific choice of the Debye 
model for the mode density which has so far been commonly used, the theory remains 
effectively unchanged as one goes over to more complicated spectrum. 

We now rewrite Eq.(12) in the form, 




(16) 
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where we use the following abbreviations, 

= x 

(17) 

u 2 = v 
and 



U\ = x 
u 2 = v 



F± = v 

F 2 = -T(X)V - V'(X) + £ eq (0 + g'(x^ ncq (t) 



(18) 



The vector u with components U\ and u 2 thus represents a point in a 2-dimensional 
'phase space' and the Eq.(16) determines the velocity at each point in this phase 
space. The conservation of points now asserts the following linear equation of motion 
for density p(u,t) in 'phase space', 



d 2 d 

— p(u, t) = - S T^F n (u, t; £ ncq , £ cq )p(w, t) , 



or more compactly 



| = -VFp. (19, 



Our next task is to find out a differential equation whose average solution is given 
by (p) where the stochastic averaging has to be performed over two noise processes £ neq 
and £ cq . To this end we note that V • F can be partitioned into two parts ; a constant 
part V • F and a fluctuating part V • Fi(t), containing these noises. Thus we write 

V • F(U, t; £neq, £eq) = V • F ( U ) + eV • F 1 ( U , t 

j Cneq; Ceq ) , (20) 

where e is a parameter (we put it as an external parameter to keep track of the order 
of the perturbation expansion in er c , where r c is the correlation time of fluctuation of 
£neq(^) ] we put e = 1 at the end of calculation) and also note that (Fi(t)) = 0. Eq.(19) 
therefore takes the following form , 

p(u,t) = (A + eAi) p(u,t) , (21) 
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where A = — V • F , A 1 = — V ■ The symbol V is used for the operator that 
differentiate everything that comes after it with respect to u. 

Making use of one of the main results for the theory of linear equation of the form 
(21) with multiplicative noise, we derive an average equation for p [(p) = p(u,t), the 
probability density of u(t) ; for details refer to Van Kampen 12 ], 

p= ^A + e 2 J™ {A^t) exp(rAo) A^t - r)> exp(-<M )} p . (22) 

The above result is based on second order cumulant expansion and is valid in the 
case that fluctuations are small but rapid and the correlation time r c is short but finite, 
i.e., 

(A^t) Ai(t')) = for \t-t'\ > t c . 

The Eq.(22) is exact in the limit correlation time r c tends to zero. Using the expressions 
for A and A 1 we obtain 

= {-V • F + e 2 / o °°rfr (V • F 1 (t) exp(-rV • F ) V • F 1 (t - r)> 

exp(rV • F )} p(u,t) . (23) 

The operator exp(— rV-F ) in the above equation provides the solution of the equation 

^ = -V. Fo f(u,t) , (24) 

(/ signifies the unperturbed part of p) which can be found explicitly in terms of char- 
acteristics curves. The equation 

u = F (u) (25) 

for fixed t determines a mapping from u{r = 0) to u(t), i.e., u — > u T with inverse 
(u T )- T = u. The solution of Eq.(24) is 



/(«,*) = /(«"*, 0) 



d{u) 



exp(-*V.F o )/(u,0) , (26) 
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d(u- r ) 



d(u) 



being a Jacobian determinant. The effect of exp(— tV ■ F ) on f(u) is as follows 

'd(u"*) 



exp(-tV-F ) /(«, 0) = /(«"*, 0) 

The above simplification when put in Eq.(23) yields 

d(-u~ 



d(u) 



(27) 



_9 
ft 



p(u,t) = V-l-F + e 2 J o 



d(u) 



{F^U^V-r-F^U-^t-T)) 

d{u) 



d(u~ T ) 



dr > p(u, t) 



(28) 



V_ r denotes differentiation with respect to U- T . We put e = 1 for the rest of the 
treatment. We now identify, 



Ui — x 
u 2 = v 

F 01 = v , F u = 
F 02 = -F(x)v - V'(x) , F 12 = + g'(x)£ Qcq (t) J 

In this notation Eq.(28) now reduces to 



(29) 







dv 



V 



d r°° 



d 



dv 



)P , (30) 



where we have used the fact that the Jacobian obey the equation 



12 



d_ 
It 



log 



d(x l , 



d(x, v) 



| B + | { -r B + ^ )} = _r 



(31) 



so that Jacobian equals to e r *. 

As a next approximation we consider the 'unpurterbed' part of Eq.(16) and take the 
variation of v during r c into account to first order in r c . Thus we have 



v + Ttv + tV'(x) 



(32) 
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Neglecting terms 0(r 2 ) Eq.(32) yields, 



dv 



- = i-rr)— + t— . 



(33) 



Taking into consideration of Eq.(33), Eq.(30) can be simplified after some algebra 
to the following form, 

^p{x,v,t) = -J^M + ^ {r(a;)u + ^'(x) - 2g\x)g"{x)I nn }p 



+ {/ ee + [(?'(a:)] 2 J nn } 



d 2 p 
dvdx 



+ {J ee - T(x)I ee + [g'(x)} 2 J nn - r(x)[^(x)] 2 7 n „ - vg\x)g"{x)I nn } 



d 2 p 
dv 2 



(34) 



(35) 



where, 

/ ee = / oo rfr(ee q (t)U(^-^))^ 
4« = dr{£ necL (t)£ necL (t - t))t 

Jee = I oo dr(Ut)Ut-r)} 

Jnn = Jo d,T{^ ne q(t)^ ne q(t t)} 

The subscripts ee and nn in the above expressions for the integrals over the correla- 
tion functions refer to equilibrium and nonequilibrium baths, respectively. In deriving 
the last Eq.(34) we have assumed that the two reservoirs are uncorrelated. Eq.(34) is 
the required generalized Fokker- Planck equation for our problem. 

In order to allow ourselves a fair comparison with Fokker-Planck equation of other 
forms 5 ' 6 ' 7 , we first turn to the diffusion terms in Eq.(34). The coefficients are coordinate 
(x) dependent. It is customary to get rid of this dependence by approximating the 
coefficients at the barrier top (say, x — 0) [one may also use mean field or steady 
state solutions of Eq.(34) obtained by neglecting the fluctuation terms and putting 
appropriate stationary condition in the diffusion coefficients]. 
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The drift term in Eq.(34) refers to the presence of a dressed potential of the form, 

R(x) = V(x) - [g'{x)f I nn 



or 

R(x) = V(x) - -7nc q ^ 2 (^) - [g'{x)\ 2 I nn . (36) 



7T 

The modification of the potential is essentially due to the nonlinear coupling of the 
system to the nonequilibrium modes. I nn is a non-Markovian small contribution and 
therefore the third term in (36) may be neglected without any loss of generality. For 
the rest of the treatment we use R(x) ~ V(x). At the vicinity of the barrier top 
x = 0, V'(x) may be approximated, as usual, by a parabolic potential, i.e., 

V(x) ^E b - \gj\x 2 (37) 



with 



V(x) ~E b - X -u? h x 2 . (38) 



For convenience, one may set g(0) = in the Taylor series expansion for g(x) (carried 
out at the barrier top x — ), without any loss of generality. And one obtains 

E b = E b (39) 

and 

^ 2 = ^ 2 + ^i[^(0)] 2 . (40) 

7T 

In the linearized description, the Fokker-Planck Eq.(34) is now reduced to the fol- 
lowing form, 

| = -,| + r P+[ r„-^]| + ^ + B ^, (4i) 

where we have used the following abbreviations; 

A = J ee - r(0)/ ee + [g\0)] 2 J nn - r(0)[</(0)] 2 / nn (42) 
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and 



B = I ee + [g'(0)} 2 I, 



(43) 



From the last two relations we have 



A = 



Jee + 9'WJnn -T(0)B 



(44) 



Defining A and B as 



one obtains 



A = ^KT and B = (3KT 



dp 
di 



dp _ 2 dp d -r^r r d 2 p 



(45) 



Identifying 



+KT 



■he + g'(0) 2 J nn 7 



d 2 p 

T(0)KT r(0)J dxdv 



n 2 = n 2 



J ee + g'(0) 2 J n 
T(0)KT 



Eq.(46) may be rewritten as, 



—p(x, v, t) = —v-^- — uj 2 x-^- + T—(vp) + ^KT 



d_ 

dt 



dx 



dv dv 



d 2 p 
dv 2 



(46) 



(47) 



+KT 



7 



ti 2 (t) 

o? r(o) 



d 2 p 
dxdv 



(48) 



Here j(t) and Q(t) are functions of time (due to the relaxation of the nonequilibrium 
modes) as defined by Eqs.(45) and (47). Or in other words nonstationary nature of the 
bath makes Cl(t) time-dependent through J nn term which is essentially a non-Markovian 
modification. 

Now the fluctuation-dissipation relations for equilibrium and nonequilibrium baths 
stated in Sec. II may be invoked. For equilibrium baths as noted earlier we have the 
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usual result; 

J ee = / dr(^ q (t)^ q (t - t)) = % q KT (49) 

J 

For the nonequilibrium version, Eq.(15) may be rearranged further to note that 

J nn = dT^ nc(i m ncq (t-r))= lncq KT(l + re-^) (50) 

J 

where r is a measure of the deviation from equilibrium at the initial instant and is 
given by r = — lj. Here U(u,t) defines the energy density distribution at 



time t. 



Using (49) and (50) we obtain from Eq.(47) 



tt 2 (t) _ r7 ncq e a* 



^ 2 ' 7c q + 7ne q [^(0)] 2 • ^ ] 

In the long time limit the relation reduces to 

Ctt^Clit) = Q . (52) 



It is interesting to note that with the replacement ~ ^% (terms are of order 

b 

1 + 0(7)) and r(0) ~ 7 one recovers the Fokker-Planck equation in the Adelman's 
form 5 (Eq.(8)). 



IV. Non-Markovian steady state Kramers rate 



We now proceed to analyze our generalized Fokker-Planck equation (48) and cal- 
culate the steady state current and the Kramers escape rate over the barrier. The 
procedure we follow in this section is similar to that of Kramers supplemented by 
Hanggi and Mojtabai's earlier analysis 6 . 

As usual we make the ansatz 

p(x, v, t) = F(x, v, t) exp 



V(x) 



KT 



(53) 
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with V(x) as approximated by a parabolic potential of the form [ see Eqs. (37-40) 

1 



with 
and 

as stated earlier. 



V(x) ~E b - ^ 2 b x 2 



*1 = ^ + 2J ^W)? 



We seek an equation for F of the form 



F(x, v, t) = F(u, t) , u = v + ax . 



(54) 



Inserting (53) and (54) in Eq.(48) we obtain 
dF f\ 1 



-,2 



(r-7)-^(r-7)^ 2 



dt 



KT 



-Axv } F 



+ [{(T- 27) -a(l + A)}v- co 2 b (l - A)x 



dF 
du 



+KT{<y + Aa) 



d 2 F 
du 2 



where 



n 2 (t) 7 



Using (51), A may be rewritten as 

A ~ 



n 2 r(o) ' 



rincqe 2* 



7e q + 7ne q [^'(0)] 2 

for r ~ 7. 

Assuming and (r — 7) to be very small we obtain 



dt du 2 



r + a(l + A) u; b 2 (l-A) ' 

V H tz : X 



T + Aa 



F + Aa 



dF 
du 



(55) 



(56) 



(57) 



(58) 
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which may be written in the form 

dF 



KT d 2 F _ dF 
dt du 2 du 



with 



a 



r + a(l + A) 



T + Aa 

and a is a solution of the quadratic equation 

a 2 (l + A) + Ta - cj 2 (1 - A) = 



(59) 



(60) 



(61) 



Since £tt^oo A = 0, the long time or steady state solution of Eq.(58) is satisfied by 

d 2 F dF 

KT ^ +a "^ = ° < 62 > 

with 



„ _ . . r + a , . 

Ctt^ooOi{t) = — = a(say) 



(63) 



Since the Eqs. (62) and (63) are identical in form to the expressions obtained in the 
usual Kramers theory one can have the usual expressions for the probability density 
p(x, v) and the current j s as 



p(x, v, oo) = N 



\ 2 rv—\a\x 

J + J dzexp 



az 



2KT 



cxp 



f + V(x) 
KT 



(64) 



with 



N 



/7rXT\2 rv-\a\x 

■ I + / az exp 



az 



2 \ 



V 2a J 'Jo r \ 2KT 

(here the subscript s in F s refers to steady state F) and 

/'+OC 3 / 2n \ 2 / Eh \ 

_ oo dvv P (^v)=N(KT)^—j exp^j 

where we have used the linearized version of V(x) near the top of the barrier at x — 0, 

V(x) =E b - l -dlx 2 , 



(65) 
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with Eb = Eb and cub is as given in Eq.(40) and N is the normalization constant. 

Employing the asymptotic distribution (just before the system is subjected to the 
shock at t — 0) of P w (x, v) for x — > — oo and at t = 0_ from p(x, v, t), where P w (x, v) = 
p(x — > —oo,v;t = 0) [see Sec. V for calculation of p{x,v,t)\, one obtains the total 
number of particles in the well, 

r +oo f +oo 2ttKT (2txKT\^ 
n a = N dv dxP w (x,v)=N . (66) 

J-oo J-oo Uo V OL J 

Here u is the frequency at the bottom of the left well. We have set the potential 
energy at the bottom of the left well equal to zero, for convenience. 

The final result for the rate of escape in the steady state is given by 

k = — = ^e~^ KT , (67) 
n a 2nujb 



where 

T\ 2 ,1* r 



A 



77) +0J 2 h 



2 



(68) 



It is evident that A is reminiscent of the 'reactive frequency' A r of Grote and Hynes 4 . 
Microscopically the non-Markovian character of the dynamics in A enters through the 
explicit structure of Y and uj b which are given by 

r = 7 e q + 7ne q [^(0)] 2 (69) 

and 

ul = ul + ^^'(0)] 2 . (70) 

71 

The appearance of the reactive frequency A is suggestive of the fact that the particle 
on the average is not moving on the bare barrier with frequency ujb but on a dressed 
barrier frequency ujb corrected by A. Pollak 8 has shown that the reactive frequency 
A is exactly an imaginary frequency of a barrier that has been modified by the bath 
modeled as a discrete set of harmonic oscillators linearly coupled to the system. The 
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effect of A is to slow down the particle by friction near the barrier. In the present 
model where the generalized Langevin equation(12) describes the motion of the particle 
over a fluctuating barrier the essential modification of A and Ub rests on the nonlinear 
coupling of the nonequilibrium relaxing modes with the system. Thus in addition to the 
properties of the bath, dynamic nature of the system-bath coupling is also significant 
in governing the barrier dynamics. We note in passing that the usual Markovian limit 
can be recovered if one puts 7 ncq = in Eq.(67) and associated quantities. 

Before closing this section one pertinent point need to be mentioned. A closer look 
into the derivation makes it clear that Eq.(67) results from an ansatz of the form 
(53) where we use V(x) in the Boltzmann factor. This choice is basically guided 
by the fact that the potential V(x) gets dressed at t = by initial excitation of 
nonequilibrium modes. This choice also makes the stationary current independent of 
position. However, if one uses the bare potential V(x) and assume a weak dependence 
of x on j s , one obtains Eq.(67) with u b in the denominator getting replaced by Uf, itself. 
The main lesson is that the modification of Kramers rate (67) is essentially due to A, 
the reactive frequency of Grote-Hynes, which has been recognized as an important 
result in view of some experimental evidence 18 of relatively weak dependence of rate 
on damping in the large friction limit. 



V. Time-dependent solution of the generalized Fokker-Planck equation ; 
nonstationary Kramers rate ; nonexponential relaxation kinetics 



We now turn to Eq.(55). Rearranging the time-dependent A-containing terms it 
may be rewritten as 



1 dF 

r¥ 



(r + a)v + 



dF KT^ F A 
du du 2 



aKT d 2 F 

r du 2 



lav 



ujIx) dF 



du 



, (71) 



where A is defined in Eq.(57). 
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Let us write 



and 



(r + a)v + ujjx 

r 

[av — <jj\x) 



= —au 



—Xu 



(72) 



(73) 



Here a is as defined in (63) and A is to be determined. 

In terms of the relations (72) and (73), Eq.(71) reduces to a more compact form. 



L?L - aU — + KT^- + A 
r dt du du 2 



aKT d 2 F ^dF' 
T du 2 du 



(74) 



Eq.(72) may be used to calculate the value of a as obtained from the solution of the 
algebraic equation 



a 2 + Ta - u 2 = 



(75) 



Only the negative root of the above equation (say a J) is the physically realizable 
solution corresponding to the steady state solution. This value of a determines uniquely 
the value of A as defined in Eq.(73) to obtain 



A = —a 



(76) 



We now seek a solution F(u,t) of Eq.(74) in the form 



F(u,t) =F s (u)e 



(77) 



where F s (u) is the steady state solution obtained in the earlier section, i.e., it satisfies 



au- 



dF), 
du 



KT 



d 2 F s 
du 2 



= . 



We require further 



(78) 



(79) 
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Substituting (77) in Eq.(74) it may be shown that the 'space' and the time part is 
separable. We obtain, 



l(ty x t _ C_ 

~Y~dt 62 ~K 



Xu 



dF s aKT d 2 F s 



du 



T du 2 



constant = -D(say) , 



where we have made use of the Eq.(78) and also 



A = Ce-i* with C = 



?~7neq 



7e q + 7ne q [^'(0)] 2 ' 



On integration over time we obtain from Eq.(80), the solution 



it 



<f)(t) = 2D-e~* 
7 



where D is determined by the initial condition. 

The time-dependent solution of Eq.(71) therefore reads as 

F(u, t) = F s (u) exp 



2DT 



-e 2 



7 



Thus the corresponding probability distribution is given by, 



p(x, v,t) — N 



{■nKT\ 2 rv-\a\x / az 

{-^r) + l dzexpi — 



exp 



KT 



e i 



2KT 



(80) 



(81) 



(82) 



(83) 



To determine D we now demand that just at the moment the system (and the 
nonthermal bath) is subjected to external excitation at t — and x — > — oo the 
distribution (75) must coincide with the usual Boltzmann distribution where the energy 
term in the Boltzmann factor in addition to usual kinetic and potential terms contains 
the initial fluctuation of energy density MA [AU = U(u, 0) — \KT\ due to excitation 
of the system at t — [see Eq.(15)]. 
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p(x, 



,v,t) — >Nl J e ~<e V 



= N 



, for(x — > — oo) . 



(84) 



The last equality demands that 



D 



j_ AU_ 
2T ~KT 



(85) 



[for the current to be coordinate independent the parabolic approximation of V(x) 
is to be used]. D is thus determined in terms of the relaxing mode parameters and 
fluctuations of the energy density distribution at t — 0. 

The time-dependent probability density therefore allows us to construct nonstation- 
ary current, 

j(t) — I dvv p(x,v,t) — j s e :r cxp ^2*) ; (86) 

J — oo 

where j s is the stationary or steady state current as derived in the last section. 
By Eq.(74) we have, 



p w (x, v) = p(x -> -oo, v,t = 0_) , 



(87) 



which was used to calculate the number of particles n a initially in the well just before 
the system was subjected to shock at t — 0. Thus non-stationary Kramers rate of 
transition is given by 



k(t) 



r 



2 



or in terms of the steady state Kramers rate k 

k{t) = k exp 



AU _ 7J 



KT 



-e 2 



(89) 
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where AU is a measure of the initial departure from the average energy density dis- 
tribution due to the preparation of the nonstationary state of the intermediate bath 
modes as a result of excitation at t — 0, and k is given by 



The above result (88) illustrates a strong nonexponential relaxation of the system 
mode undergoing a nonstationary activated rate process. The origin of this is an initial 
preparation of nonequilibrium mode density distribution (with a deviation AU) which 
eventually relaxes to an equilibrium distribution. Eq. (88) implies that the initial 
transient rate is different from the asymptotic steady state Kramers rate. What is 
immediately apparent is that the sign of AU[— U(uj,0) — \KT] determines whether 
the initial rate will be faster or slower than the steady state rate. When AU is negative, 
i.e., the contribution of thermal energy dominates, the initial rate of thermal activation 
of the reaction co-ordinate gets enhanced as a consequence. On the other hand, when 
the sudden excitation of the nonequilibrium modes provides a positive deviation AU, 
the initial rate of activation becomes slower. This is because there likely to exist some 
time lag for the nonthermal energy gained by the few nonequilibrium modes by sudden 
excitation to be distributed over a range before it become available to the reaction 
co-ordinate as thermal energy for activation. 

It is also interesting to consider the zero and high temperature limits. When 
T — > both the steady state Kramers rate k as well as the time-dependent factor 
exp[— -^e _7t//2 ] goes to zero. If T = 0, then k(t) is zero at all time. However, it seems 
intuitively that there should be a transient period during which the rate is finite. It 
may be noted that since the relaxation of the nonequilibrium bath modes (following 
the sudden excitation) is very slow compared to the rate of activation process, the par- 
ticle undergoing barrier crossing cannot 'sense' this transient (ideally if the relaxation 
to equilibrium is adiabatic, i.e., the thermalization of the initial departure AU is very 




(90) 
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slow, there should be no transient). We believe that the distinct separation of the two 
time-scales implied in the dynamics makes the transient unobservable. An interplay of 
overlapping time-scales pertaining to the relaxation of the bath and the activation of 
the system may give rise to transients in k(t) at T = 0. Evidently this is outside the 
scope of the present treatment. 

When T is very high such that \KT far exceeds U(u, 0) the initial rate gets strongly 
enhanced (since AU is negative) and the time- dependent exponential factor becomes 
roughly independent of temperature. In the limit t — > oo or AU — > we recover steady 
state Kramers rate, as expected. 

The activation rate is thus consequently modified which effectively incorporates a 
secondary relaxation kinetics. The quasi-thermal excitations decay on the time scale 

which is well separated from other internal time scales of the thermal bath. The 
dynamic nature of the coupling between the system and the nonequilibrium modes is 
responsible for fluctuating barrier. A closer look into the origin of the non-exponential 
kinetics makes it clear that the spiritual root of D-term is essentially the A-containing 
term in Eq.(71) or term in Eq.(46) which is a non-Markovian contribution. We 
thus identify the non-exponential relaxation of the system mode as a typical non- 
Markovian dynamical feature. In the case of very small 7 one naturally recovers the 
exponential relaxation and Arrhenious rate of activation of the usual kinetic scheme. 

A relevant pertinent point regarding some of the related works need be considered 
here. Generalized Langevin equation (GLE) has been widely employed in various con- 
texts, e.g., in the description of reactions in liquids. A search for realistic models began 
with the realization that friction exerted by the solvent on the solute is space depen- 
dent. A formally consistent approach to the problem of space and time dependent 
friction had been introduced early by Lindenberg and co-workers 20 ' 21 . Carmeli and 
Nitzan 22 have also derived a stochastic dynamical equation which is a generalization 
of GLE to the case of space and time dependent friction. Pollak and Berezhkorskii 25 
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have demonstrated that the space and time-dependent friction model is identical to a 
multidimensional anisotropic but Markovian friction problem in which the reaction co- 
ordinate is coupled to an additional co-ordinate which is governed by a Langevin type 
equation. A theory for treating spatially dependent friction in the classical activated 
rate processes has been considered and following the method of Pollak an effective 
Grote-Hynes reactive frequency for this case has been obtained as a transcendental 
equation 23 . More recently a general theory for thermally activated rate constants in- 
fluenced by spatially dependent and time correlated friction 24 has been proposed. 

While in the above problems one is concerned with the space and time dependent 
friction, which is essentially a characteristic of the solvent mode structure, in the present 
problem we deal with effect of a secondary relaxation of intermediate oscillator modes 
(following an initial excitation) on the primary kinetics of the system mode. The mode 
density function due to initial excitation differs from its equilibrium value - a feature 
which is marked in the nonequilibrium fluctuation-dissipation relation. Thus the expo- 
nential relaxation in Eq.(81) is not be confused with the exponential time- dependent 
friction employed in earlier instances. The origin of these two exponential terms are 
fundamentally different. The non-exponential kinetics is essentially an offshoot of a 
dynamic modification of the fluctuation-dissipation theorem appropriately carried over 
to a nonstationary regime. This nonequilibrium nature of activated process is reflected 
in the nonstationary kinetics that we derive here. 

The non - exponential relaxation kinetics had been explored earlier in different 
occasions in relation to disordered systems 13 , viscous liquids 19 , oxygen binding to 
haemoglobin 16 , where phenomeno logical fluctuating barrier models have been employed 
(barriers arising from the collective motions of many degrees of freedom). The present 
model although oversimplified in many respects captures the essential nature of influ- 
ence of an initial non-thermal mode density distribution on the relaxation kinetics of 
the system. 
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VI. Conclusions 



In conclusion, we consider a simple microscopic system-nonequilibrium bath model 
to simulate nonstationary thermally activated processes. The nonequilibrium bath is 
effectively realized in terms of a semi-infinite dimensional broad-band reservoir which 
is subsequently kept in contact with a thermal reservoir which allows the nonthermal 
bath to relax with a characteristic time. A systematic separation of timescales is then 
used to construct the appropriate Langevin equation for the particle, which is nonlinear 
and non-Markovian in character. Based on a strategy of Van Kampen's expansion in 
er c of the relevant physical quantity where e is the strength and r c is the correlation 
time of fluctuations of the relaxing modes, we show that this Langevin equation can 
be recast into the form of a generalized Fokker-Planck equation, when the correlation 
time is short but finite. Adelman's form of the Fokker-Planck equation [ Eq.(8) ] as 
well as the standard Markovian description can be recovered in the appropriate limits. 
We now summarize the main conclusions of this study: 

(i) The model proposed here captures the essential features of Langevin dynamics 
with a fluctuating barrier. The present approach is equipped to deal with situations 
both in the non-stationary short time as well as stationary long time regimes. The 
origin of the short time non-exponential kinetics can be traced back in a non- stationary 
fluctuation-dissipation theorem. 

(ii) We derive the expression for the steady state Kramers escape rate in the non- 
Markovian case and show that the Grote-Hynes 'reactive frequency' can be realized 
explicitly in terms of the microscopic parameters of the nonequilibrium relaxing modes 
and their arbitrary dynamic coupling to the system mode. 

(iii) The central result of this paper is the derivation of a nonstationary Kramers rate 
in closed analytic form. This essentially illustrates the influence of an initial excitation 
and subsequent relaxation of the nonequilibrium bath modes on the system degree of 
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freedom undergoing an activated process. The system mode is shown to follow strong 
non-exponential kinetics. 

The model considered in the present paper may be realized in a guest-host system 
embedded in a lattice where the immediate local neighborhood of the guest comprises 
intermediate oscillator modes whereas the lattice plays the role of a thermal bath. 
Appropriately identified reaction co-ordinate coupled to other degrees of freedom in 
a molecule embedded in a matrix may be another worthwhile candidate for such a 
scheme. 

Although simple, the model thus allows us explicit solutions and in view of the 
prototypical role played by the present model in several earlier investigations, we hope 
that the conclusions drawn here will find applications in some related experiments of 
physics and chemistry of complex systems. 
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